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We present a detailed assessment of how well a future Electron-Ion Collider could constrain helicity 
parton distributions in the nucleon and, therefore, unveil the role of the intrinsic spin of quarks and 
gluons in the proton’s spin budget. Any remaining deficit in this decomposition will provide the 
best indirect constraint on the contribution due to the total orbital angular momenta of quarks 
and gluons. Specifically, all our studies are performed in the context of global QCD analyses based 
on realistic pseudo-data and in the light of the most recent data obtained from polarized proton- 
proton collisions at BNL-RHIC that have provided evidence for a significant gluon polarization in 
the accessible, albeit limited range of momentum fractions. We also present projections on what can 
be achieved on the gluon’s helicity distribution by the end of BNL-RHIC operations. All estimates 
of current and projected uncertainties are performed with the robust Lagrange multiplier technique. 
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I. INTRODUCTION AND MOTIVATION 

The exploration of the nucleon’s inner structure and 
the interactions among its constituents at high energies 
has been the main protagonist in the quest for a quan¬ 
titative picture of matter at the most elementary level 
during the last fifty years. In spite of an impressive 
wealth of achievements and discoveries, many fundamen¬ 
tal questions still remain unanswered today. The origin 
of the proton spin is a remarkable example for one of 
these compelling questions still driving the field of Nu¬ 
clear Physics. Inspired by the Quark Parton Model, it 
was originally thought to arise solely from the intrinsic 
spins of the proton’s valence up and down quarks, but 
this naive view was denied by results from polarized deep 
inelastic scattering (DIS) experiments in the late eight¬ 
ies [H and remained under experimental and theoretical 
scrutiny and debate ever since. 

The pi'OTOsed Electron-Ion Collider (EIC) project in 
the U.S. ^ versatile machine designed to explore 

nuclei and polarized light ions at the energy and inten¬ 
sity frontier with the precision of electromagnetic and 
electroweak probes, will be the next milestone in the 
quest for gaining a deeper insight into the nucleon’s in¬ 
ner workings. Among its remarkable and unprecedented 
capabilities, an EIC will be able to disclose definitively 
the role of quark and gluon spins in the proton down to 
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fractions a; of a few hundred thousandths of the nucleon’s 
momentum and in a wide range of resolution scales set 
by virtuality Q of the probing photon in inclusive deep- 
inelastic scattering (DIS). Once the intrinsic spin contri¬ 
butions are being tightly constrained, an EIC will deliver 
early on a first quantitative, though indirect, estimate on 
the combined orbital angular momenta (0AM) of quarks 
and gluons needed in the decomposition of the proton 
spin. Of course, an EIC is also aiming at directly access¬ 
ing 0AM [2^ which is, however, experimentally much 
more challenging than inclusive DIS and, hence, will re¬ 
quire some years of running and additional theoretical 
groundwork to bear fruit. 

In the phenomenologically very successful framework 
of perturbative QCD (pQCD) information on how mo¬ 
mentum and spin of quarks and gluons are apportioned in 
the nucleon factorizes in a universal, process-independent 
way from calculable, short-distance partonic scattering 
cross sections Q. In case the incident beams are both 
longitudinally polarized in the scattering process, helic¬ 
ity parton distributions (PDFs) Af{x, Q^) are accessible, 
where / denotes the different (anti-)quark flavors q (q) 
or the gluon g. For factorization to be a viable approach, 
the external scale Q characterizing the scattering needs 
to be large enough, say, above about 1 — 2 GeV. Upon 
integration over all momentum fractions x, 

A/(Q2)= ['dxAf{x,Q^) (I) 

Jo 

the helicity PDFs contain the desired information enter- 
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ing in the decomposition of the proton’s spin 

2 = 2 ^ ( 2 ) 

where £((5^) = + Lg{Q'^) is the 

total contribution from 0AM. The sum of quark and an¬ 
tiquark densities in Eq. ([5]) can be combined in the quark 
singlet contribution AE((5^). It should be noted that 
the apparently simple relation ([^l was subject to quite 
some debate and controversy in recent years. The main 
problems are that the decomposition ([2) is not unique 
and that the separation of the gluonic total angular mo¬ 
mentum into Ag{Q^) and Lg{Q^) was not thought to 
be possible in a gauge invariant manner. It is now un¬ 
derstood 0 that a particular choice of decomposition is 
essentially a matter of taste as long as each component 
can be determined, at least in principle, experimentally 
or from lattice QCD calculations. Each variant covers 
complementary aspects of a complex bound-state system 
such as the nucleon. 

In this paper, we are primarily interested in providing 
the best possible quantitative assessment of the impact 
a future EIC would have on determinations of the quark 
singlet and gluon helicity densities and their contribu¬ 
tions to the proton spin in Eq. ©■ The theoretical frame¬ 
work for our survey is based on an earlier study Q where 
we have performed a series of global QCD analyses of he¬ 
licity PDFs at next-to-leading order (NLO) accuracy in¬ 
cluding sets of realistic mock EIC data for both polarized 
DIS and semi-inclusive DIS (SIDIS). These data were 
generated and properly randomized within one sigma un¬ 
certainties using state-of-art spin-dependent PDFs 00 
based on experimental information from polarized fixed 
target DIS and SIDIS and proton-proton collisions at 
BNL-RHIC available at that time. 

Since then, many new and important results by the 
BNL-RHIC experiments [TTI - [l^ have been reported, that 
have changed substantially our perception of the gluon 
helicity distribution El-' QCD analyses of first RHIC 
spin results i in were dominated by single-inclusive 
pion production data, obtained at rather low transverse 
momentum scales pr, and suggested very little or no 
gluon polarization in the range 0.05 ^ x ^ 0.2 predom¬ 
inantly probed by data, albeit within very large uncer¬ 
tainties. In contrast, the most recent and very precise jet 
production data n from RHIC, taken at larger scales 
Pt than for pions, provided for the first time clear evi¬ 
dence for a sizable gluon polarization [13; very similar 
results have been reported in an independent global QCD 
analysis in Ref. [l^. It is important to notice that both 
sets of recent inclusive measurements at RHIC which are 
highly sensitive to the gluon helicity PDF, jets [l^ and 
pions [l^, are fully consistent with each other if one 
properly takes into account the different scales pt and 
ranges of momentum fractions x probed by the respec¬ 
tive data. This is automatically guaranteed and achieved 
in the framework of a global QCD analysis at NLO ac¬ 


curacy based on exact, unabridged NLO expressions for 
the underlying spin-dependent pp cross sections Eng. 

In light of these new results from RHIC pertaining to 
the relevance of gluons in the nucleon’s spin decomposi¬ 
tion it is of critical importance to reexamine our previous 
impact study 0 to update the case for a future polarized 
DIS program at an EIC. Apart from including the latest 
RHIC spin data we will also refine and expand our study 
in many other important aspects. To establish a mean¬ 
ingful baseline for estimating the impact of future EIC 
DIS data, we first present projections of what will likely 
be the final word from the RHIC spin program with re¬ 
spect to Ag{x, Q"^) and its integral Ag{Q^). To this end, 
we add realistic projections for yet to be published in¬ 
clusive jet and pion data from polarized pp collisions at 
200 and 510 GeV [2l| to the global QCD analysis frame¬ 
work of Ref. [13. First, preliminary results for these 
measurements are fully consistent with the most recent 
global fits of helicity PDFs [HO- Next, we update and 
include the sets of mock polarized DIS data to reflect the 
latest collision energies envisioned for the eRHIC option 
of an EIC Q. Throughout, we pay special attention to 
properly estimate the impact of the different data sets in 
the various regions of parton momentum fraction x. 

Once we have completed our extensive studies of the 
prospects for the gluon helicity density Ag{x, Q^) and 
its integral Ag{Q^) from both the RHIC spin program 
and a future EIC, we will consider the quark spin con¬ 
tribution AE(a;,(5^) summed over all flavors. Finally, we 
will examine what can be gleaned from these two projec¬ 
tions for the total parton 0AM contribution C{Q^) when 
combined with the spin decomposition @. 

As a technical improvement, we introduce in all our 
impact studies a more robust treatment of uncertain¬ 
ties. Specifically, we now estimate uncertainties with the 
Lagrange multiplier technique [l^ considering dynami¬ 
cal tolerances corresponding to the 90% confidence level 
(C.L.) [13 as it is custo mary in most current fits of unpo¬ 
larized parton densities |24l |. rather than adopting some 
ad hoc fixed Ay;^ criterion as in our previous analyses 
[M3. This new procedure takes into account more pre¬ 
cisely the role of the different observables in constraining 
parton densities at different kinematic regions in x and 


We have to postpone a new detailed study of the fla¬ 
vor separated helicity quark distributions Af{x, Q^) as a 
proper global QCD analysis of projected polarized SIDIS 
data with identified pions and kaons for an EIC critically 
depends on having available a reliable set of parton-to- 
hadron fragmentation functions (EEs) and their uncer¬ 
tainties. At this point, only pion EEs have been updated 
recently 2^ and kaon EEs, important for determining 
the strangeness polarization, are still work in progress 
[ 13 . Also, RHIC spin data, this time for the longitudinal 
single spin asymmetry of bosons [^ are expected to 
play a significant role in establishing a meaningful base¬ 
line for impact studies of flavor separated Af{x, from 
SIDIS measurements at an EIC. However, the final data 
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sets from RHIC are not yet available. Finally, we note 
that for this study we focus in general on the so far un¬ 
measured region of small momentum fractions and, cor¬ 
respondingly, low-to-medium values of at an EIC. The 
prospects of DIS measurements with an EIC in the re gion 
where Q ~ Mw was studied in some detail in Ref. 
based on NLO expressions for charged current DIS given, 
for instance, in (28l| . 

In the remainder of the paper, we first briefly discuss 
the sets of pseudo-data generated to study the impact 
of both polarized DIS at an EIC and upcoming RHIC 
pp data in determinations of the gluon and quark sin¬ 
glet helicity distribution and their respective integrals. 
In Secs. Ill and IV we present the results of our global 
QCD analyses at NLO accuracy based on including these 
sets of data one-by-one. First, we focus on what can be 
expected for Ag{x, Q^) by the end of the RHIC spin pro¬ 
gram and, next, we include also the mock EIC data. Our 
main results are best presented in terms of “running x- 
integrals” for Ag (in Sec. HI), AE, and the difference 
between 1/2 and ^AS-I-Ag (in Sec. IV), that represents 
roughly how much room is left for 0AM to saturate the 
spin sum rule. Finally, we summarize our results and 
present our conclusions. 


II. SIMULATED DATA FOR POLARIZED DIS 
AT AN EIC 

Inclusive DIS with longitudinally polarized leptons and 
nucleons can be expressed by the structure function 
gi{x,Q'^). For energy scales Q well below the mass of 
the electroweak gauge bosons and Z, the scattering 
is mediated by the exchange of a virtual photon and both 
kinematic variables x and can be straightforwardly re¬ 
constructed from measuring the energy and angle of the 
deflected lepton. 

Sensitivity to the gluon helicity distribution Ag{x, Q^) 
in DIS is manifold: directly, through higher order QCD 
corrections to gi , which are sub-leading compared to the 
dominant quark contribution and, hence, small and diffi¬ 
cult to utilize in an analysis, or, indirectly, through scal¬ 
ing violations, that is variations of the structure function 
gi{x, Q^) with scale Q for fixed values of x. The smaller x 
the more pronounced are the scaling violations, which are 
usually expressed as the logarithmic derivative of the DIS 
structure function, i.e., dgi{x,Q^)/d It should 

be noted that corresponding results for unpolarized DIS 
from the DESY-HERA experiments [1^ are utilized in 
all current extractions of helicity-averaged PDFs and pro¬ 
vide by far the best constraint on the gluon density below 
about X ~ 0.05 [l^. Through evolution, Ag{x, Q^) is also 
correlated with the quark singlet helicity distribution, at 
small X mainly to the sea component, which itself is very 
poorly constrained below the range of existing polarized 
DIS data covering x > 0.004 for > 1 GeV^ [ll,US]- As 
we shall see below, this also leads to a significant uncer¬ 
tainty in current estimates of gi{x, Q^) and, hence, in the 


quark singlet AE(Q^) entering ([^l, which only an EIC 
can finally resolve. At variance with unpolarized PDFs, 
where the scale evolution at small x is driven by gluons, 
quark and gluon helicity PDFs are equally important as 
they exhibit a similar, less singular a: —>■ 0 behavior in 
all evolution kernels . The different evolution of he¬ 
licity and helicity-averaged distributions also leads to a 
strong suppression of gluon and quark polarizations as 
a: —^ 0 and, consequently, of all spin asymmetries in that 
kinematic region. 

Thus, in order to get an accurate picture of the gluon 
and quark helicity densities and, in particular, their x- 
integrals m at an EIC it is not only crucial to have good 
precision in a wide range of parton momentum fractions x 
to reduce extrapolation uncertainties from the currently 
unmeasured x region in DIS, but also to cover for any 
given X the largest possible range in the photon virtuality 
to determine Ag from scaling violations. The require¬ 
ment of being able to reach x values well below of what 
has been achieved so far 0,113 for virtualities safely in 
the DIS regime, say, above at least 1 GeV^, clearly calls 
for the largest conceivable center-of-mass system (c.m.s.) 
energies y/s at a future ep collider in order to make any 
significant impact. Furthermore, at small x the reach in 
to study scaling violations, which are only logarithmic 
in energy, is kinematically limited to = sxy where y 
is the fractional energy of the virtual photon, y is con¬ 
strained by the increasing depolarization of the virtual 
photon from below, which should not exceed about 90%, 
and, from above, by the energy of the scattered lepton, 
which should be not too low to be still reliably measur¬ 
able in the detector. In practice, O.OI < g < 0.95 appears 
to be a reasonable choice as the DESY HERA collider ex¬ 
periments even used the range 0.005 ^ y ^ 0.95 for parts 
of their physics program. 

TABLE I: Combinations of electron and proton energies of 
the current eRHIC design [J as used in our analyses, the 
corresponding c.m.s. energies, and the range in x accessible 
for = 1 GeV^ and assuming 0.01 < y < 0.95. 


Ee X Ep 

[GeV] 

[GeV] 

^min 

for y — 0.95 

^max 

for y — 0.01 

15 X 100 

77.5 

1.8 X 10”'‘ 

1.7 X lO""" 

15 X 250 

122.7 

7.0 X 10"® 

6.7 X 10“® 

20 X 250 

141.4 

5.3 X 10"® 

5.0 X 10“® 


In our previous study Q, we adopted c.m.s. energies 
from ^/s ~ 45 GeV to 71 GeV corresponding to an elec¬ 
tron beam energy of 5 GeV for the BNL eRHIC design of 
an initial version of an EIC at that time. As a potential 
later upgrade we also considered = 20 GeV leading 
to a maximum ^/s = 141.4 GeV. The latest eRHIC de¬ 
sign evolved considerably and now has a default electron 
beam energy of about 15 GeV but is also capable to run 
with both lower and higher lepton energies at similar lu¬ 
minosities 01 from day one. When combined with the 
existing RHIC proton beam of up to Ep = 250 GeV this 
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set-up not only extends the reach towards small x by al¬ 
most two decades compared to what is known today but 
also makes the detection of the scattered electron at large 
y much easier than for a 5 GeV electron beam. Table |T] 
summarizes the three different c.m.s. energies we are go¬ 
ing to consider in this study along with the accessible 
range in x, [xmin, a^max], for = 1 GeV^, assuming the 
same standard DIS cuts as in our previous analysis Q: 

> lGeV^ > 10GeV^ and 0.01 <y< 0.95. As 
can be seen, the improvements in the kinematic range 
for both X and from i/s = 122.7 to 141.4 GeV corre¬ 
sponding to Ee = 15 GeV and 20 GeV, respectively, are 
not too significant. 

In order to generate the mock polarized DIS data for 
the energies shown in Tab. HI we proceed very much in 
the same way as was detailed in Ref. Q. The statistical 
accuracy for each of the generated data sets is scaled 
to correspond to a very modest accumulated integrated 
luminosity of 10 fb”^, equivalent to at most a few months 
of operations for the anticipated luminosities for eRHlG 
0 and 50% efficiency in the data taking. We assume 
70% polarization for both the electron and the proton 
beam. Each of the three new sets adds about 70 points 
to the existing suite of fixed target DIS spin asymmetry 
data, distributed logarithmically in 4 [5] bins per decade 
in [a;]. As in Ref. @ , the actual data used in our global 
analyses below are not the generated ones but theoretical 
estimates of the spin asymmetry at NLO accuracy based 
on the latest DSSV helicity densities [13, reflecting the 
same relative statistical accuracy in each x, bin as 
the Monte Carlo data, and having their central values 
randomized within one-sigma uncertainties. 

In any case, the DIS measurements at an EIC will 
be mainly limited by systematical errors due to detector 
performance, beam polarization and luminosity measure¬ 
ments, and the unfolding of QED radiative corrections, 
which we all assume to be controlled to a percent level 
accuracy. This will be necessary for achieving meaning¬ 
ful constraints on the helicity PDFs as we have already 
demonstrated in our previous analysis, see, e.g.. Fig. 6 in 
Ref. We note that the typical size of the DIS spin 
asymmetries at the lowest values of x can be as small as a 
few times 10“^. Such measurements are already routinely 
performed at RHIC in case of double spin asymmetries 
for single-inclusive pion production at transverse mo¬ 
menta pt of a few GeV. 

To illustrate the main idea of the planned measure¬ 
ments in polarized DIS, the enhanced kinematic coverage 
at an EIC, and the current uncertainties which prevent 
us from having a clear picture of how quarks and glu¬ 
ons generate the proton’s spin, we start with presenting 
the scale dependence of the gluon helicity density. Fig¬ 
ure [T] shows the shape of Ag{x,Q'^) as extracted from 
the latest global QCD analysis to all presently available 
data in DIS, SIDIS, and pp collisions [13] (henceforth de¬ 
noted as “DSSV 2014 best fit”). The plot illustrates how 
the QCD scale evolution quickly broadens the a;-shape of 
Ag(x,Q^) with increasing Q^, which is initially largely 



FIG. 1: [color online] QCD scale evolntion of the optimum 
DSSV 2014 gluon helicity distribution as a function of 
the momentum fraction x for various values of (dotted 
curves). The dashed lines delineate the minimal x accessible 
at an EIC in the x — plane for various c.m.s. energies 
and y — 0.95. The line for ^/s = 44.5 GeV corresponds to 
an energy considered in our previous paper utilizing only 
a 5 GeV electron beam at eRHIC. The shaded area labeled 
“RHIC” shows approximately the region where Ap(x, Q^) is 
constrained best by published RHIC pp data [l^ . [l^ . 


concentrated at rather high momentum fractions a; at a 
scale of = 1 GeV^, and moves its peak rapidly towards 
smaller values of x. The main constraint on Ag{x, Q^) 
so far is derived from the latest RHIC inclusive jet and 
7 r° [H, [l^ which span, as is indicated in Fig. [U only a 
rather limited range in both x and the latter be¬ 
ing approximately determined by the measured range 
of PT (a few GeV to about 20-30 GeV). Clearly, noth¬ 
ing is known about the low momentum fraction tail of 
Ag{x, Q^) beyond the part of it generated radiatively by 
evolution from lower and larger x. In the next Sec¬ 
tion, we will elaborate in detail on the present uncertain¬ 
ties for Ag(x,Q^), which appear to very substantial, in 
particular, in the unmeasured small x region. In general, 
spin asymmetries sensitive to small x are very difficult 
to seize experimentally as the unpolarized gluon density 
rises sharply as a; —I 0, leading to a more and more sup¬ 
pressed gluon polarization Ag/g, likewise for sea quark 
polarizations. 

To emphasize the importance of running an EIC at the 
largest possible c.m.s. energy, we also overlay in Fig. [TJ 
the kinematic lower limits in x and for various val¬ 
ues of ^/s and assuming y < 0.95, see Tab. H] The line 
for ^/s = 44.5 GeV corresponds to the no longer consid¬ 
ered lowest c.m.s. energy of our previous study Q; as 
was mentioned above, the new design for eRHIC allows 
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FIG. 2: [color online] Projected EIC data for the structure 
function gi{x,Q^) for the different combinations of electron 
and proton energies in Tab. H] Constants are added to gi to 
separate the different x bins and multiple data points in the 
same {x, Q^) bin are slightly displaced horizontally. The solid 
lines are obtained for the optimum DSSV fit of 2014 0 and 
the shaded bands illustrate the 90% C.L. uncertainties due to 
variations in the gluon helicity density. The shaded region in 
the lower left corner illustrates the {x,Q^) region covered by 
present fixed target data. 


to cover significantly lower values of x from the very be¬ 
ginning of operations. 

Figure [2] illustrates our updated simulated data sets 
for inclusive polarized DIS at an EIC for the three differ¬ 
ent choices of c.m.s. energies listed in Tab. ID The solid 
lines reflect the expectations from the best fit of DSSV 
2014 [0 by extrapolating their results outside the exper¬ 
imentally constrained x and range. The shaded bands 
illustrate the uncertainty estimates correspon ding to the 
90% C.L. variations of lS.g{x, given in Ref. [I 3 , which 
cover a very significant spread below about x ~ 0.01; 
see also Fig. 1 in Ref. [oj. The error bars for the EIC 
pseudo-data were determined as described above and in 
Ref. Q and reflect the expected statistical accuracy for a 
modest integrated luminosity of lOfb”^, 70% beam po¬ 
larization, and 50% efficiency in the data taking. We re¬ 
call that all currently available polarized DIS data cover 
only the lower left corner in Fig. [5] with the smallest a;, 



FIG. 3: [color online] The polarized DIS structure function 
g\{x,Q^) at — lOGeV^ as a function of x computed with 
the optimum DSSV 2014 helicity PDFs [itI ] (solid line). The 
dotted curves represent alternative fits within 90% G.L. un¬ 
certainties. The dashed and dot-dashed lines show the effects 
of the scale evolution from — 1 GeV^ to 100 GeV^. The 
points illustrate typical uncertainties and the kinematic reach 
of projected EIC data for the three different c.m.s. energies 
listed in Tab. H] 

X ~ 3.6 X 10“^, being reached by the recent COMPASS 
data [ 3 ^ for ~ 1 GeV^. As can be seen, in the kine¬ 
matic region already covered well by present fixed target 
data, X > 0.01, the remaining uncertainties in gi{x,Q‘^) 
are very small. For smaller x, the precision of the pro¬ 
jected EIC data is significantly better than current un¬ 
certainties and these measurements will be the decisive 
factor in future global fits as we shall illustrate in the 
next Section. 

One notices the rather modest scaling violations 
d 5 i(x, Q^)/dlnQ^ for the optimum DSSV 2014 fit 
throughout the entire x and range shown in Fig. [2l 
in particular, if compared to similar plots for the unpo¬ 
larized DIS structure functions [ 2 ^. On the one hand, 
this is due to the less singular scale evolution for helicity 
PDFs at small x, and, on the other hand, there is also 
a potential delicate cancellation with the quark helicity 
PDFs, which, as Ag itself, are not bound to be positive 
definite and, in addition, can have different signs for dif¬ 
ferent flavors. Therefore, alternative fits, like those for 
Ag shown in Fig. [1] will all exhibit somewhat different 
patterns of scaling violations than the optimum DSSV 
2014 fit. 

As we shall see next, our current ignorance of the small 
X behavior of helicity quark densities also imposes a sig¬ 
nificant uncertainty on expectations for gi{x,Q^) in the 
EIC regime. In Fig. [3] we present the DIS structure func- 
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tion gi{x,Q‘^) (solid line) and 90% C.L. estimates of its 
uncertainties (dotted lines) as a function of the momen¬ 
tum fraction a; at = lOGeV^. Unlike in Fig. O the 
alternative fits at 90% C.L. now include combined vari¬ 
ations of quark and gluo n helicity PDFs away from the 
DSSV 2014 best fit [l7l | which lead to uncertainties at 
least twice as large as for the variations just based on Ag 
shown in Fig. [5] We note that throughout our paper the 
allowed ranges of variations at 90% C.L. are determined 
for each of the shown results by the robust Lagrange 
multiplier technique and dynamic tolerances for the ap¬ 
propriate increase in the of fh® similar to what is 
done in most of the recent PDF fits 2^ . 

To illustrate once again the accuracy of future mea¬ 
surements at an EIC, we also show here a few repre¬ 
sentative projected data points taken from Fig. [5] in the 
relevant regime around 10 CeV^ for the three differ¬ 
ent c.m.s. energies we consider. Clearly, measurements 
of gi{x,Q^) at small x will dramatically reduce the un¬ 
certainties in the quark helicity PDFs and, indirectly, 
through the coupled QCD scale evolution of quarks and 
gluons also in Ag{x,Q‘^). At any given x, scaling viola¬ 
tions for gi{x,Q^) will further constrain Ag{x,Q^). As 
was already shown in Fig. [2l they are numerically not 
very pronounced for the optimum DSSV 2014 fit, which 
can be also inferred from Fig. [31 where we show gi{x, Q^) 
at = 1 and 100 CeV^ in addition to our default scale 
of 10 CeV^. However, each of the alternative fits exhibits 
a somewhat different dependence driven by the uncer¬ 
tainties in the x shapes of the quark and gluon densities. 
For X > 0.01, the scale dependence of gi{x,Q^) in the 
range from = 1 to 100 CeV^ is typically larger than 
the uncertainty on gi{x,Q'^) from present data. 


III. PRESENT STATUS OF Ag AND IMPACT 
OF PROJECTED RHIC AND EIC DATA 

Before addressing the question of how precisely an EIC 
will constrain the total gluon and quark polarizations in 
the spin decomposition ([3]), and, indirectly, also the total 
0AM £, it is important to first make a precise assessment 
of how well these quantities are expected to be known by 
the end of the current experimental programs, in partic¬ 
ular, RHIC spin. This will set the best possible baseline 
to judge on the impact a future EIC could have in the 
field of QCD spin physics. 

Different indicators and measures can be adopted to 
quantify how well the gluon helicity density and the re¬ 
sulting contribution Ag{Q'^) to the proton’s spin are con¬ 
strained by data. The standard way to study uncertain¬ 
ties as a function of the parton’s momentum fraction x 
at a given in a global QCD fit to all available data 
is certainly the most obvious possibility, however, it nei¬ 
ther provides an immediate idea of the accuracy for the 
phenomenologically interesting x-integral © that is the 
focus of our study, nor does it indicate the relevance of 
the different regions in x probed by the different experi- 



FIG. 4: [color online] The running integral of the gluon he¬ 
licity distribution at = 10 GeV^ as a function of Xmin ac¬ 
cording to the DSSV 2014 analysis 0 (solid line) and 90% 
C.L. uncertainty estimates (shaded bands) based on global 
QCD fits utilizing different sets of existing and projected pp 
and EIC data (see text). 


ments used in the fit. 

Instead, we choose to present most of our results in 
terms of the “running integral” of, for instance, the gluon 
helicity density, defined analogously to Eq. o as 

Ag{Q'^,Xmin) = f dxAg{x,Q'^), (3) 

^min 

which represents the share of the proton spin ([3]) from 
gluons as a function of the lower integration limit Xmin- 
Its uncertainty takes into account the non-trivial corre¬ 
lations between the different regions of x contributing to 
(O . By varying Xmin in © , one can explore how low in x 
- or, alternatively, how high in -^5 - one likely needs to go 
with future experiments to reduce x —> 0 extrapolation 
uncertainties to a level small enough to make meaningful 
statements about how gluons and quarks in the proton 
make up its spin. To study the important question of 
the convergence of © with Xmin in more detail, we will 
also compute the contributions to © from different bins 
[xinini a^max] in X in case of Ag. 

To estimate the impact of past, current, and future 
data sets on Ag and AE we proceed in steps. To this 
end, we will present uncertainty estimates for various 
running integrals by including different data sets one- 
by-one into our global analysis framework. As we have 
mentioned already, to demonstrate the impact an EIC 
will have on Ag in the future, we should take into account 
the experimental information that is expected to become 
available soon from the RHIC spin program. Essentially, 
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the RHIC running focusing on longitudinally polarized 
double-spin asymmetries has concluded in 2015, and sev¬ 
eral high impact data sets are forthcoming, some of which 
have been presented in preliminary form at conferences 
recently [^. Since the expected statistical accuracy of 
these measurements is already known, we can proceed 
by generating mock RHIC pp data with the proper un¬ 
certainties to estimate their impact on helicity PDFs, in 
particular, the gluon. 

Since extending the reach towards smaller values of x 
is the most important asset to arrive at a more solid esti¬ 
mate for the integral Ag(Q^), we only focus on upcoming 
RHIC measurements which have the best potential to do 
so, i.e., inclusive jet and neutral pion data from STAR 
and PHENIX, respectively, taken at the highest c.m.s. 
energy of 510 GeV for pp collisions at RHIC. In particu¬ 
lar, spin asymmetries at forward rapidities will allow one 
to probe x values down to about a few times 10“^, signifi¬ 
cantly below the currently covered range in x indicated in 
Fig.[TJ More specifically, our projections comprise double 
spin asymmetries for 

• inclusive jet data at mid rapidity, 

• inclusive data at mid rapidity, and 

• inclusive 7r° data at forward rapidities . 

In Fig. m we show the running integral ([3]) of the gluon 
helicity distribution at = lOGeV^ down to Xmin = 
10“®. Concentrating first on the optimum DSSV 2014 
fit (solid line), one notices that the integral saturates 
and reaches more than 90% of its total value Ag(Q^ = 
lOGeV^) already for Xmm — 10“^, suggesting that most 
of the gluon spin contribution to the sum ([2]) stems from 
momentum fractions above 10“^. This is in line with 
common expectations that Ag{x,Q^) behaves like 
xg{x, Q^) at sufficiently small values of x. The preferred 
value for the gluon spin contribution at = 10 GeV^ in 
the DSSV 2014 analysis [a turns out to be rather large, 
about 0.37, corresponding to roughly 70% of the proton 
spin. However, neither the preferred small x behavior, 
nor the large value for Ag{Q^ = lOGeV^), are based on 
actual experimental constraints in the QCD analysis as 
is illustrated by the various uncertainty bands shown in 
Fig. m Here and in the following, all uncertainty bands 
correspond to 90% G.L. estimates. 

The outermost light shaded bands in Fig. |4] represent 
a global QCD fit based solely on existing polarized fixed 
target DIS and SIDIS data [ij. One can safely conclude 
that these data do not constrain the gluon helicity density 
or its running integral by any means. The importance of 
the current RH IC p olarized pp data for inclusive jet and 
7 r° production [lll - llq in constraining Ag{x, Q"^) down 
to about X ~ O.OI is exemplified by the second level of 
shaded bands. Still, for smaller values of Xmin possible 
variations in the truncated integral away from its best fit 
value quickly become very large due to the flexible func¬ 
tional form adopted in the DSSV framework. At least, 


the analysis suggests at 90% C.L., that the truncated in¬ 
tegral ([3|) stays positive down to a;min — 2 x 10“^, but 
nothing can be concluded about the full integral Ag{Q^), 
not even its sign. These uncertainties correspond to those 
obtained in the DSSV 2014 global analysis [13 and rep¬ 
resent our current knowledge. 

Next, we include the projections for the upcoming 
RHIC pp data as listed above into the global analysis 
framework to estimate what can be expected to be known 
about Ag{Q^) by the end of RHIC operations and be¬ 
fore an EIC will turn on. Any impact of a future EIC 
needs to be judged with respect to this baseline. As can 
be seen, RHIC can be expected to further constrain the 
gluon helicity distribution and, in particular, its running 
integral, mainly by ruling out very extreme, positive or 
negative, variations at small momentum fractions. This 
will be mainly achieved by measurements at forward ra¬ 
pidities. Even though these data will have no sensitivity 
below X ~ 10“^, they considerably limit possible varia¬ 
tions within the given very flexible functional form of the 
DSSV analysis. Together with the decreasing weight of 
the region of very small x in the running integral ([3|) , the 
error on the integral down to Xmin — I0“® is significantly 
reduced but it still covers about twice the proton spin of 
1/2 and ranges roughly from —0.25 to 0.75. 

We are now in a position to focus on the impact of the 
generated EIC data shown in Fig.jH To this end, we suc¬ 
cessively add the three sets of mock data given in Tab. U 
starting from the lowest c.m.s. energy of ^/s = 77.5 GeV 
to the global set of data, including the RHIC projections 
just discussed. The series of the innermost dark shaded 
bands in Fig. |4] represents the corresponding uncertainty 
estimates. The results clearly show that for momentum 
fractions x > 0.05, DIS, SIDIS, and RHIC data, as in¬ 
cluded in the DSSV 2014 analysis [13 , constrain the run¬ 
ning integral adequately well whereas an EIC adds very 
little new here, apart from a crucial, independent check of 
what has been learnt so far in that particular kinematic 
region. 

The constraining power of an EIC starts to become in¬ 
creasingly noticeable below x ~ 0.01. Already at around 
ccmin ~ 10“^, EIC DIS data are expected to reduce the 
uncertainties by approximately a factor of 8 as compared 
to the DSSV 2014 estimate, and by a factor of 4 relative 
to the projections based on future RHIC data. Below 
Xmin — 10“^, EIC data will constrain the flexible func¬ 
tional form for Ag{x, adopted both in the DSSV and 
our analyses in such a way that uncertainties in the run¬ 
ning integral remain constant, excluding extreme varia¬ 
tions in the gluon helicity density at small x still allowed 
in fits based solely on current or projected RHIC data. 

Using only EIC data with the lowest c.m.s. energy con¬ 
sidered in our analysis, would lead to a determination of 
the gluon spin contribution to (I2|) to within a relative 
uncertainty of about 25%. Including data at higher y/s 
would further reduce relative uncertainties in Ag{Q^) to 
a level of about 10%. We stress again, that for all our 
studies, the EIC mock data are generated and random- 
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FIG. 5: [color online] Contributions to Ag{Q^) in ([2]) at = 
10 GeV^ from different x intervals (see text). The solid and 
dashed lines rraresent the results from the DSSV 2014 0 
and 2008 [^, [l0| best fits, respectively, which mainly differ 
by latest RHIC pp data [la . [l^ . The various shaded bands 
reflect 90% C.L. uncertainty estimates for both fits and the 
impact of forthcoming RHIC data. The dotted line indicates 
the cumulative contribution to Ag{Q^) for DSSV 2014. 


ized around the small x extrapolation of the DSSV 2014 
best fit (solid line in Fig. 2]). An EIC would be able to 
either verify if gluons indeed contribute only very little 
to the running integral below Xmin — 10“^, as is usually 
expected (3^ . or to determine decisively if they exhibit 
any different behavior at small x. In the latter case, de¬ 
pending on at which aJmin the running integral starts to 
saturate to its full value, DIS data at the largest possible 
c.m.s. energy are potentially most relevant to have. 

The relevance that different regions in x have in the 
running integral shown above is analyzed more closely in 
Figs. |S] and m Figure|S]illustrates the comparative impact 
of gluons to the spin sum ([^l in different bins of x^ four 
per decade, and the corresponding uncertainties down to 
10“^ expected to be accessible with forthcoming RHIC 
pp data. The solid (dashed) lines in each bin are the 
estimates computed with the DSSV 2014 [l^ and 2008 
0 fits, respectively. The dotted line gives the cumulative 
contribution to the running integral m as was shown in 
Fig- HI The arrow indicates the value of full integral 
= lOGeV^). In addition, 90% C.L. uncertainty 
estimates are shown for each individual bin. We notice 
that the partial contributions of each bin in x to the 
running integral ([3]) and the corresponding variations are 
very strongly correlated to those of the neighboring bins. 
As a consequence, the sum of the individual uncertainties 
from each bin can exceed the uncertainty of the running 
integral shown in Fig. 01 Nevertheless, they correctly 


FIG. 6: [color online] As in Fig. [5] but now focusing on the 
small X extrapolation and the impact of the projected polar¬ 
ized DIS data sets shown in Fig.[2]on the 90% C.L. uncertainty 
estimates. 


reflect the relative impact of the different bins in x. In 
each bin the uncertainties reflect the synergy of all the 
different data sets included in the global analysis. 

The outermost shaded bands in Fig. [5] correspond to 
the DSSV 2008 analysis, the middle ones to the most re¬ 
cent 2014 fit. The optimum gluon helicity distribution 
in the 2008 fit had an almost vanishing contribution to 
the proton spin sum rule @ and was characterized by 
a node in the x region probed by available RHIC data 
at that time. These features can be easily gathered also 
from Fig.jS] The current 2014 best fit based on the latest 
RHIC pp data has a much larger and positive Ag{x, Q^), 
cf. Fig. dJ As can be inferred from Fig. [SJ only the x 
region probed by the pp data contributes significantly to 
the integral. As is expected from Fig. 01 uncertainties for 
the old DSSV 2008 analysis are significantly larger than 
those for the 2014 fit in all bins of x, but both fits can¬ 
not constrain contributions to the running integral be¬ 
low X ^ 0.01 in any meaningful way. Here, uncertainties 
sharply increase with decreasing x. The innermost, dark 
shaded uncertainty bands demonstrate the constraining 
power of the projected RHIC data which now, unlike for 
the previous results, exhibit decreasing uncertainties with 
smaller values of x. 

In Fig. we extend the x range down to 10“® and 
include bin-by-bin uncertainty estimates based on the 
three sets of projected DIS data at an EIC. Compared to 
Fig. 01 we now use only 2 bins per decade in x. Below 
X ^ 10“^ the contributions to the integral in each bin 
remain small for the DSSV 2014 best fit (solid line), fluc¬ 
tuating around zero, but the corresponding uncertainties 
continue to increase for both the 2008 and 2014 fits. Un- 
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certainty estimates including also the projected RHIC 
data only slightly increase towards smaller x reflecting 
the improvements in the running integral already ob¬ 
served in Fig. m The constraining power of the EIC DIS 
data is very significant in the entire x range and appears 
to be roughly a constant factor of 2 — 3 with respect to 
the bin-by-bin estimates including the projected RHIC 
data. This is somewhat less than what was obtained for 
the running integral at low ccmin, suggesting that scal¬ 
ing violations are more powerful when considered over 
an extended range of momentum fractions rather than 
in a small bin in x. In the latter case much fewer data 
points are actually providing a strong constraint on al¬ 
lowed variations in any given bin. Another important 
factor to consider are the mentioned sizable bin-by-bin 
correlations. Similarly, the impact of including more DIS 
sets at different c.m.s. energies is reduced for the bin-by¬ 
bin studies as compared to running integral shown in 
Fig.H 


IV. STATUS AND PROSPECTS FOR AS AND 
THE TOTAL OAM CONTRIBUTION 

Similarly to what has been discussed in the previous 
Section in connection with the gluon helicity density, the 
running integral of the quark singlet AE(a;, repre¬ 
sents the intrinsic spin contribution of all quark flavors 
in the decomposition of the proton spin ([5]). Thanks to 
the direct coupling of the quarks to the probing virtual 
photon in DIS, AE(x, is much better constrained by 
present fixed target data than the gluon helicity distri¬ 
bution which only enters indirectly through QCD scale 
evolution or as an 0{as) correction. Since AE(a;, 
and Ag{x, are coupled through the singlet evolution 
equations, any constraint from data on either of the two 
distributions impacts also the other one. 

An extraction of the quark singlet from DIS data 
on gi{x,Q^) also requires to determine simultaneously 
two additional flavor non-singlet distributions, which, if 
needed, can be all recast into the total contributions 
from u, d, and s quarks, i.e., Au + Au, Ad + Ad, and 
As -|- As (here we ignore for simplicity any contribution 
from charm and bottom quarks, which play no role for 
all currently available data). 

The x-integrals of the two non-singlet combinations are 
usually assumed to be related to the hyperon decay con¬ 
stants F and D within some uncertainties, which provide 
some indirect constraint for the currently unmeasured 
small X region below a few times 10“^. Since we only 
wish to focus on the quark flavor singlet in this paper, 
we adopt these constraints in the same way as was done 
in the DSSV global analyses. A more thorough analysis 
of flavor separated quark helicity densities and their un¬ 
certainties will be conducted once a fully updated suite 
of reliable sets of fragmentation functions becomes avail¬ 
able to revisit our previous impact study of SIDIS data 
in Ref. Q. Then we will also explore how well an EIC 



FIG. 7: [color online] Similar to Fig.[4]but now for the running 
integral of the quark singlet helicity density AE. The solid 
line corresponds to the optimum fit of DSSV 2014 [T^. 90% 
C.L. uncertainty estimates (shaded bands) are shown for the 
DSSV 2008 and 2014 analyses and after including the different 
sets of projected EIC data. 


can challenge the constraints imposed on the quark sec¬ 
tor by the hyperon decays, one of which being related 
to the Bjorken sum rule, the other, more important one 
mainly to the amount of strangeness polarization in the 
nucleon. 

In Fig. [3 we show 90% C.L. estimates for the running 
integral of AE(x, Q^) as a function of Xmin for = 
10 GeV^ for fits including different sets of existing and 
projected data. One should notice that the vertical axis 
here covers only half of the range shown for the running 
integral of Ag in Fig. |4l Also, AE((5^) enters the spin 
sum rule ([2]) with a factor of 1/2 relative to the gluon 
spin contribution. 

The outermost shaded band represents uncertainties as 
present in the original DSSV global analysis from 2008 
Q- They appear to be very significant for Xmin ^ 10“^. 
As usual, the solid line shows the optimum fit of DSSV 
2014 extrapolated down in x. The corresponding uncer¬ 
tainties are much reduced as compared to the 2008 anal¬ 
ysis due to including additional DIS data from the COM¬ 
PASS collaboration and, indirectly, through constraints 
on Ag from RHIC pp data [l^ IT^ . 

The set of three innermost shaded bands illustrates 
the significant impact of an EIC from a series of global 
fits that successively include the projected DIS data 
sets starting from the one corresponding to the lowest 
c.m.s. energy. The addition of the sets with increasing 
c.m.s. energy has less impact on the uncertainties than 
for the gluon helicity distribution shown in Fig. 21 For 
a^min = 10“® and with an EIC, one expects from our stud- 
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FIG. 8: [color online] As in Fig. [7|but now showing the contri¬ 
bution from the combined quark and gluon 0AM C{Q^) us¬ 
ing Eq. © and the results obtained for the running integrals 
of Ag and AE. As before, 90% C.L. uncertainty estimates 
(shaded bands) illustrate the impact of different projected 
EIC data sets. 

ies to control the value of AE to within about 15% which 
is somewhat worse than what can be achieved for Ag. 
Also, the convergence of the running integral for AE for 
the DSSV 2014 best fit is much slower than for the gluon 
density shown in Fig. 31 This can be mainly attributed 
to the small x behavior of the strangeness helicity dis¬ 
tribution and the above mentioned imposed constraints 
from the hyperon decay constants F and D. An EIC 
will be able to verify the validity of these assumptions 
by determining the sea quark distributions, in particu¬ 
lar, As{x,Q‘^), down to unprecedentedly small values of 
X from a series of SIDIS measurements with identified 
pions and kaons. Once this is achieved, we expect that a 
future combined global analysis of SIDIS and DIS data 
will also lead to a much improved uncertainty estimate on 
the integral of AE compared to what is presented here. 

As we have demonstrated in Figs. [Hand [71 an EIC will 
deliver a precise picture of the intrinsic spin contributions 
of quarks and gluons to within at least 10 — 15% rel¬ 
ative uncertainties. This information can be used along 
with the proton’s spin decomposition (©) to estimate how 
much is left for the combined contribution C{Q'^) from 
the orbital motion of quarks and gluons. This is illus¬ 
trated in terms of the running integral of C{Q‘^) in Fig. [51 
along with our usual set of uncertainty estimates based 
on the different sets of existing and projected experimen¬ 
tal data. 

The solid line gives the estimate for C{Q^) for = 
lOGeV^ based on extrapolating the best fit results of 
DSSV 2014 for Ag(x,Q^) and AT,{x,Q^) down to Smin 
and subtracting them off the total proton spin of 1/2. 


As it turns at, at this particular value of Q^, the DSSV 
2014 result converges to basically a zero net contribution 
from 0AM to the spin sum rule (©) but within huge un¬ 
certainties. The latter is given for both the DSSV 2008 
and 2014 global fits by the outermost and middle shaded 
bands, respectively. As one can already anticipate from 
our studies performed in Figs. [Hand[71 an EIC will yield 
an excellent indirect constraint on C{Q‘^) by combining 
the then available precise information on Ag and AE. 
As can be inferred from the corresponding uncertainty 
estimates (set of innermost shaded bands), one can ex¬ 
pect to constrain the net 0AM contribution of quarks 
and gluons to within about ±0.05 when all three sets of 
projected DIS data are combined. Of course, only a set 
of dedicated measurements at an EIC can reveal a de¬ 
tailed, hopefully flavor separated, quantitative picture of 
the orbital motion of quarks and gluons. For instance, 
even if the net 0AM contribution turns out to be small 
as for the DSSV 2014 best fit, there might be significant 
cancellations among the individual quark flavors and the 
gluon. In addition, Ag{Q^) in (©) evolves logarithmically 
with such that any increase has to be compensated 
by a corresponding decrease of £((5^), and vice versa, to 
always arrive at 1/2. We recall that AE(Q^) is a renor¬ 
malization group invariant at LO and, in the MS scheme, 
evolves only very slowly with Q^. However, separating 
quark and gluon 0AM in Eq. (©) experimentally involves 
determinations of twist-three generalized parton as well 
as quantum phase-space Wigner distributions which will 
be very challenging and still needs further theoretical 
work and perhaps lattice QCD studies, see, for instance. 
Refs. 0 and [s^. 


V. SUMMARY AND CONCLUSIONS 

The Electron Ion Collider project constitutes a versa¬ 
tile and vast program to considerably deepen our knowl¬ 
edge of the inner workings of nucleons and nuclei and the 
underlying dynamics and interactions of quarks and glu¬ 
ons as governed by Quantum Chromodynamics. We have 
presented a detailed account of what can be achieved in 
the field of Spin Physics, in particular, the complex inter¬ 
play of quark and gluon spins and their orbital motions 
yielding the known spin 1/2 of the proton, for which we 
still lack a detailed quantitative understanding despite 
more than thirty years of intense research. 

In the present study, we have considerably updated 
and refined our previous analysis on the impact of future 
EIC polarized DIS data in connection to our understand¬ 
ing of the gluon spin and its share in the proton spin 
budget. Throughout, we have explored the consequences 
of the much larger degree of gluon polarization in the nu¬ 
cleon than it has hitherto been acknowledged, as implied 
by the most recent RHIC results. In addition, we have 
used realistic projections for several forthcoming sets of 
RHIC data to estimate to what extent they further con¬ 
strain the gluon helicity distribution by the end of the 
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RHIC program. This baseline was used to quantify the 
expected impact of a future polarized DIS program at an 
EIC. 

To this end, we have focused on providing accurate 
estimates of uncertainties for the running integrals as 
a function of the minimum momentum fraction for the 
gluon and quark singlet helicity distributions. Results 
were obtained in the context of global QCD analyses 
based on current and various projected sets of data and 
using the robust Lagrange multiplier technique for all er¬ 
ror estimations. In case of the gluon helicity distribution 
we also provided detailed uncertainty studies bin-by-bin 
in momentum fraction to establish where the running 
integral receives major contributions and to determine 
below which value of x one can expect the integral to 
converge. 

We find that the RHIC spin program still has quite 
some potential in reducing current uncertainties in the 
gluon helicity distribution by ruling out extreme varia¬ 
tions in its shape. Nevertheless, one cannot expect to ar¬ 
rive at any meaningful and reliable estimate of the gluon 
contribution in the proton spin decomposition due to re¬ 
maining substantial extrapolation uncertainties from the 
unmeasured region of small momentum fractions. Even 
under very conservative assumptions for its performance, 
an EIC with a sufficiently high enough center-of-mass 
system energy to reach deep into the so far unexplored 
kinematic regime, will decisively constrain both the gluon 


and the singlet quark helicity densities and, at last, ele¬ 
vate the field of QCD spin physics into an precision era. 
For the corresponding a;-integrals relevant for the pro¬ 
ton’s spin decomposition we expect to achieve relative 
uncertainties of about 10%. For the first time, this will 
also lead to a precise indirect estimate of the remaining 
contribution to the proton spin from the combined orbital 
motion of quarks and gluons, well before an EIC will be 
able to determine them directly. Estimates obtained by 
extrapolating our current best knowledge of helicity den¬ 
sities, without EIC data, to small momentum fractions 
point to a surprisingly small net orbital angular momen¬ 
tum in the proton’s spin budget at intermediate scales 
Q^, albeit within very large uncertainties. 
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